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Properties of Vector Preisach Models 
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Abstract — This paper discusses rotational anisotropy and 
rotational accommodation of magnetic particle tape. These effects 
have a performance impact during the reading and writing of the 
recording process. We introduce the reduced vector model as the 
basis for the computations. Rotational magnetization models must 
accurately compute the anisotropic characteristics of ellipsoidally 
magnetizable media. An ellipticity factor is derived for these 
media that computes the two-dimensional magnetization 
trajectory for all applied fields. An orientation correction must be 
applied to the computed rotational magnetization. For isotropic 
materials, an orientation correction has been developed and 
presented. For anisotropic materials, an orientation correction is 
introduced. 

Index Terms — Preisach modeling, reduced model, rotational 
anisotropy, rotational accommodation, rotational orientation. 

I. INTRODUCTION 

T HIS paper presents three extensions of the vector Preisach 
model: the reduced vector Preisach model, RVPM, which 
eliminates the need for a rotational correction in the 
SVPM, and extension of the orientation correction to 
anisotropic media, and a technique for handling vector 
accommodation. 

The accurate characterization of magnetic processes 
requires a vector model of magnetic hysteresis. The classical 
Preisach model cannot adequately represent vector magnetic 
processes since it is inherently a scalar model. Several authors 
have modified the scalar Preisach model to include the vector 
features of a magnetic medium [1-4]. The Simplified Vector 
Preisach Model (SVPM) [1] was developed for computing the 
vector magnetization in response to a vector-applied field. 

The SVPM is a coupled-hysteron model that exhibits both 
the saturation property and the loss property [1], The vector 
magnetization is computed from the integration of the product 
of a state vector and a Preisach function and then performing 
the rotational correction. The state vector is computed using 
the selection rules determined by the applied field. The 
rotational correction term provides the cross-axis coupling 
effect of applied fields. This cross-axis coupling term makes it 
impossible to obtain an analytical closed form solution for the 
magnetic susceptibility. 

The Reduced Vector Preisach Model (RVPM) does not 
require the rotational correction. The developed vector model 
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uses modified selection rules for state vector calculation. 

Using this new model as a basis to characterize the medium, 
measurements and computations of the angle of magnetization 
for a magnetic particle (MP) tape rotated in an applied field 
have been made. This anisotropic material is ellipsoidally 
magnetizeable, that is, when subject to a constant magnitude 
rotating field has a magnetization trajectory that is nearly 
elliptical [11]. This model has an orientation error which has 
been corrected for isotropic media. We introduce a correction 
for anisotropic media. 

An ellipticity factor has been computed for ellipsoidally 
anisotropic media that accurately computes the two- 
dimensional magnetization for large applied fields (50% 
greater than the hard axis coercivity of the material). A 
procedure has been derived and is presented for computing the 
magnetization trajectories for small applied fields. The 
measured magnetization trajectories (polar plots of the 
magnetization magnitude and angle) of elliptically 
magnetizable material changes shape from ellipses for small- 
applied fields to circles for large applied fields. The measured 
ellipticity factor, X, for rotational applied fields approaches 
one asymptotically as the applied field increases and the 
magnetization trajectories become circles. 

The phenomenon of accommodation has been presented and 
discussed for scalar Preisach models [9] and a scalar 
differential equation model has been introduced [10] to 
compute accommodation. Rotational magnetization 
measurements have demonstrated vector accommodation for 
ellipsoidally magnetizable media [11]. From the magnetization 
trajectories, accommodation is demonstrated for 360° of 
rotation. The radial change in a magnetization trajectory from 
the start to the end of applied field rotation is a measure of the 
accommodation. As the rotating applied field increases, the 
magnetization trajectories start with shapes of ellipses and 
exponentially approach the shape of a circle. 

n. Simplified Vector Preisach Model 

The SVPM computes the normalized irreversible 
magnetization components as the product of the rotational 
correction R(l x , I y , I z ), and the basic Preisach integrals 

m ,j = R(l x JyJ, )l j , for j = x, y, or z, (1) 

The output from the basic Preisach integrals Ij are computed as 
I j = j| Qj p{uj , Vj )dUjdVj , for j = x, y, or z, (2) 

v i <u i 

where the up and down switching fields are m ; and v ; , the 
normalized Preisach function is p, and the state function is Qj. 
The state function Qj is determined by selection rules. These 
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selection rules are summarized in Table I for the 2-D model 
and are generalized for the 3-D model in Table IQ. 

The rotational correction is given by 



It was shown [1] that for any set of values of /,, the rotation 
correction ensures the correct magnitude of the magnetization 

and 1 < R < -fi . 

The normalized irreversible magnetization is computed as 
the vector sum of three basic Preisach models as 

mj = RI . (4) 

The irreversible magnetization can be computed as the vector 
sum 

M, =M s S«m 1 , (5) 

where M s is the saturation magnetization and S is the material 
squareness matrix defined as 

S x 0 O' 

s= 0 S y 0 . (6) 

0 0 S z 

Thus, to simulate anisotropic media, this model can allow 
different values for the S’s along each of the axes and 
parameters in the basic Preisach models. 

The normalized reversible magnetization components can 
be computed as 

m Ri = a j+fiPj )- «/-/(- h j ) * foT J = * y- °r 4 (7) 
where the a }± variables can be implemented with either a state- 
independent, magnetization-dependent, or state-dependent 
reversible magnetization as in the case of the scalar models. 
The non-linear function /(// ; ) is defined as 

f(pj)=l-e-*>, ( 8 ) 

where £ is a model parameter that needs to be identified. 

The reversible magnetization can be computed as a vector 
sum, 

M r =(l-S)M s m R , (9) 

where S is the squareness matrix defined in (6). 

Then the total magnetization can be expressed as 

M t =M,+M r . (10) 

m. Reduced Vector Preisach Model 

The Reduced Vector Preisach Model computes the 
normalized irreversible magnetization components using the 
basic Preisach integrals 

m Ij = ttQRjPfaj’VjfajdVj . for; = x, y, or z, (11) 

v i <u i 

where the up and down switching fields are and v the 
normalized Preisach function is p, and the state function is 

Qr /. 

The state vector, Q Rl , for the RVPM differs from Qj for the 
SVPM, and is computed by new selection rules summarized in 
Table I for the 2-D case. These rules differ from the SVPM 
case only at the comers. The subscript d is used for the 


direction in which Q d is being computed. The subscript c is 
used to indicate the cross direction. The new selection rules 
are defined such that no rotational correction is required for 
computing the magnetization [5]. 

Similarly, the irreversible magnetization, the reversible 
magnetization component and the total magnetization are 
computed using (5), (9) and (10) respectively. 

The application of selection rules shown in Table II shows 
that at any point on the Preisach hyperplane, the sum of the 
square of the Cartesian components of the state vector obeys 

Q 2 +Q 2 +Q 2 =l. ( 12 ) 

For the ellipsoidal magnetization behavior, the major 
remanence path must satisfy 

ml+ml+ml= 1. (13) 

For sufficiently strong fields, the normalized reversible 
magnetization satisfies 

m 2 Rx +m 2 Ry +m 2 Rz = 1. (14) 

Therefore, the normalized total magnetization can be 
expressed as 

mr x + m j y + m} z =1. (15) 

The important properties of the RVPM can be summarized 
as: 

The model is applicable to anisotropic media as well as 
isotropic media. 

The model reduces to the scalar model when the 
applied field and the initial magnetization lie along only 
one of the principal axes, whereby, the vector model 
will have all the properties of the scalar moving model. 
In the presence of large fields, the normalized irreversible 
magnetization and the reversible magnetization each trace out 
ellipses whose major axes are at right angle to each other for 
both isotropic and anisotropic media. 

IV. Reduced Vector Preisach Model Simulations 

The RVPM is applied to an isotropic magnetic medium. For 
an isotropic medium, values for O f and O k are equal, and 

negligible compared to the average critical field H k . The 
material parameters used for simulations are: O ( . = O k =165, 

H k = 633,5= 0.57, M s = 0.014127,0= 33332.81, £ = 
0.0009. 

Simulations were carried out for the orthogonal component 
H r of the applied field of 1.1ft* , 1.4ft* and 4 ft* . The vector 
DOK model [7] is implemented based on the RVPM using the 
cobweb method [6] for computation speed. Figure 1 shows the 
magnetization angle versus the applied field angle. It can be 
seen that the magnetization ratchets as the applied field rotates. 
For very large applied fields, the magnetization angle becomes 
equal to the applied field angle for all values. 

Figure 2 shows the locus of the magnetization as the applied 
field is rotated. It is seen that as the field increases, the curves 
become rounder. The flattening of the loci close to 135° and 
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315° is a discretization error caused by the jump in the 
magnetization angle for the respective applied fields. 

For isotropic media, since the magnetization rotates faster 
than the applied field, the rotations of the magnetization need 
to be corrected using the correction rules defined in [8]. 
Applying these corrections for the applied fields as in Fig. 1, 
the corrected magnetization angle vs. applied field angle is 
plotted as shown in Fig. 3. Figure 4 shows the magnetization 
loci with corrections for an applied field rotation of 360°. It is 
seen that the curves are more rounded. 



Fig. 1. Magnetization angle vs. applied field angle for applied fields 
of 1.1 A* (solid line), 1.4 (dotted line), and 4 (dash-dot line). 



Fig. 2. A plot of the corrected magnetization angle vs. applied field 
angle for the same conditions as in Fig. 1 . 



Fig. 3. Locus of magnetization for the same set of applied field 
rotations as shown in Fig. 1. 



Fig. 4. A plot of the corrected magnetization angle vs. applied field 
angle for the same conditions as in Fig. 1. 

V. Orientation 

A simplified vector model was introduced [3] that is the 
basis for the isotropic model. Two orthogonal Preisach planes, 
each with its own Preisach function, are defined. For the easy 
axis Preisach function, seven Preisach parameters (saturation 
magnetization, squareness, average critical field, zero field 
susceptibility, standard deviation of interaction field, standard 
deviation of critical field, and moving parameter) are 
identified for the isotropic material. The Preisach parameters 
used for each orthogonal Preisach plane are set equal to this 
set of identified parameters, thereby, making the model 
isotropic. The model then computes the irreversible 
magnetization incorporating a rotational correction factor. 
From the irreversible magnetization, the reversible 
magnetization is computed. The irreversible and reversible 
magnetizations are vectorally added together to obtain the total 
two-dimensional magnetization. This isotropic computation 
contains an orientation error as shown in Fig. 1. For isotropic 
material, the magnetization should follow the applied field 
after the magnetization begins to rotate. Applying the 
orientation correction [2] removes this error as shown in Fig 2. 
This paper extends the correction from the isotropic model to 
the anisotropic model. 

The two-dimensional vector magnetization of MP tape is 
computed using a cobweb representation of a two-dimensional 
Gaussian Preisach function [12] and [6], where the first 
dimension is the X-Preisach plane and the second dimension is 
the Y-Preisach plane. Figure 5 shows a cobweb for the X- 
Preisach plane. Outside the outer ring of the cobweb, the 
Preisach function is considered to have a negligible value. 
This Preisach function is defined in each plane using seven 
Preisach parameters identified for the MP tape in each plane. 
The process of computing the magnetization for the 
anisotropic MP tape consists of first computing the isotropic 
magnetization, applying an isotropic orientation correction to 
this magnetization, and then applying an anisotropic correction 
to obtain the anisotropic magnetization. 
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Fig. 5. Cobweb Preisach x-plane centered at the coercivity of a 
material, H*. Outer ring of cobweb is assumed to be at least 3ai 
and 3 at from H^. 



Fig. 6. . MPT computed anisotropic magnetization-applied field angle 
curves computed from the corrected isotropic magnetization-applied field 
angle curves. The dotted line is a reference line where the magnetization 
angle equals the applied field angle. 


To compute the isotropic magnetization for the MP tape, 
the Preisach parameters for the hard axis are set equal to the 
parameters identified for the easy axis. Using this isotropic 
model, the isotropic magnetization is computed. Application 
of the orientation correction removes the orientation error as 
shown in Fig. 2. 

The magnetization-applied field angle curves measured for 
anisotropic MP tape are double periodic [11] An anisotropic 
correction factor must be applied to the model results to obtain 
the double-periodic oscillating characteristic of anisotropic 
media. Since this oscillating characteristic is sinusoidal in 
nature, a sinusoidal correction may be derived. Let 9 h be the 
angle of applied field after applying the orientation correction. 
Then 


X = sin(fy) 


(16) 


and 

Y = cos(0 A ) (17) 


which matches the measured curve, except for applied fields 
between plus and minus the coercivity, H c . Figures 6 and 7 
show the computed, using (18) and (19), and the measured 
anisotropic magnetization-applied field angle curves, 
respectively. 

The measured magnetization-applied field curves are 
computed from measured x- and y-component magnetization 
for the applied fields p^H = 0.5, 0.25, and 1.0 T. The curves 
are double periodic. The angle of magnetization exceeds the 
angle of applied field before 180° and before 360“ (0“), which 
correspond to the negative and positive easy-axis directions of 
the MP tape. This result happens, because the angle of 
magnetization lies between the easy axis and the angle of 
applied field. When the x-component of the applied field 
exceeds the coercivity, H c , the magnetization tends toward the 
closest easy axis, resulting 


are the sinusoidal representations of the anisotropic applied 
field, X, and magnetization, Y, respectively. There may be a 
required linear shift associated with (16) and (17). Let 0 H and 
0m be the anisotropic angles of applied field and 
magnetization, respectively. Then, 


} H -X +A 0 H , 

(18) 

[ay' 

atan2 +A0^f, 

\ X J 

(19) 


where A d H and A0 M are the required associated linear shifts, 
and X is the anisotropic correction factor. The variables, A 0 H 
and A 0 M , respectively, determine the amount of horizontal shift 
of 9n and the amount of vertical shift of 0 M . The amplitude of 
the sinusoidal curve, Y ', in (19) is adjusted by X and generates a 
double periodic curve between 0° and 360°, 



Fig. 7. MPT measured anisotropic magnetization-applied field angle 
curves for three applied fields. The dotted line is a reference line where 
the magnetization angle equals the applied field angle. 
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in a possible angle of magnetization greater than the angle of 
applied field. As the applied field increases, the curves tend to 
straighten out and lie along the dashed 45° line, where the 
angles of magnetization and applied field are equal. 

Figure 8 shows the anisotropic correction, X, used in (19) to 
adjust the amplitude of the sinusoidal to fit the measured 
curves in Fig. 7. This curve is a smooth ascending curve to 
H c - Above H c , this curve is a smooth descending curve 
approaching the value of one asymptotically as the applied 
field increases. The flat appearance of this curve in the 
neighborhood of He is due to the discretation size of the 
applied field, goH = 0.05 T. 



Kg. 8. Measured ellipticity factor for MP tape over the applied field 
range // = 0 T to 1.3 T. 

VI. Orientation Extension 

The SVPM and the RVPM each produce a magnetization 
angle-applied field angle curve with a four-fold periodicity as 
shown in Fig. 1. Since the angle of magnetization actually 
follows the angle of applied field uniformly, an orientation 
correction factor is applied by the model to ensure that this 
characteristic is correctly produced. The orientation correction 
factor removes the four-fold periodicity from 0° to 360°, but it 
also removes the two-fold anisotropic periodicity that is 
characteristic of anisotropic media shown in Fig. 7. 

As the applied field increases from zero, the magnetization 
starts and remains at zero until the applied field begins to 
sweep through the media’s hard-axis Preisach function. At 
this point, the magnetization begins to rotate and follow the 
applied field. As the applied field continues to rotate, the 
magnetization continues to rotate as more hysterons switch. 

Because the media are anisotropic, the hysterons along the 
easy axis of the media switch easier than the hysterons along 
the hard axis switch, thereby introducing a sinusoidal 
periodicity into the magnetization. The magnetization always 
lies between the applied field and the easy axis. Rotating the 
applied field from 0° to 360°, the magnetization follows the 
applied field until the easy-axis component of the 
magnetization switches from positive to negative. When the 
easy-axis component of the magnetization becomes negative. 


the magnetization lies between the applied field and the 
negative easy axis (180°); the magnetization now leads the 
applied field as shown in Fig. 7. The magnetization leads the 
applied field until the easy-axis component of the applied field 
becomes negative, and the magnetization returns to follow the 
applied field. 

To implement a two-dimensional rotational anisotropic 
model, one must consider how much of the hard axis Preisach 
function is swept by the applied field as shown in Fig. 9 for a 
positive applied field. In Fig. 9, the black area of the Preisach 
function is the area swept by the applied field. When the 
applied field is less than the Preisach function, none of the 
Preisach function is swept, and the Preisach function in Fig. 9 
is all white. When the applied field is greater than the 
Preisach function, the entire Preisach function in Fig. 9 is 
black. The areas of the Preisach function that are black 
correspond to magnetic particle whose magnetization rotates, 
whereas, the white areas correspond to magnetization that 



Kg. 9. Preisach function along hard axis of MP tape with applied field, 
ytoH = 0.15 T. Magnetization rotates in black area; magnetization 
oscillates in white area as the applied field rotates. 

oscillates. Because the hard axis applied field varies from 
positive to negative during rotation of 360°, there is a s imil ar 
sweeping of the hard axis Preisach function for negative 
applied fields. 

The computation of anisotropic magnetization follows the 
process depicted by the flow diagram shown in Fig. 10. When 
the applied field does not sweep the Preisach function, P R = 0, 
there is no rotation of the magnetization, and the magnetization 
is computed by the RVPM without an ellipticity correction. 
When the applied field completely sweeps the Preisach 
function, Pr = 1, there is complete rotation of the 
magnetization, and an ellipticity correction is applied to the 
entire RVPM computed magnetization. When the Preisach 
function is partially swept by the applied field, 0< P R <1, a 
portion of the magnetization is rotating, C R , and a portion of 
the magnetization is oscillating, C 0 - The 
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Hr = Applied Field 



Hg. 10. Flow chart to compute rotating and oscillating magnetization. 

oscillating magnetization is computed by multiplying the 
RVPM magnetization by the coefficient, Co- The rotating 
magnetization is computed by multiplying the RVPM 
magnetization by the ellipticity correction factor and then by 
the coefficient, C R . The total magnetization is then determined 
by vectorally adding the oscillation magnetization and the 
rotational magnetization as shown in Fig. 11. 

VII. Ellipticity Factor 

The ellipticity factor is the ratio of the change in 
magnetization along the easy axis to the change in 
magnetization along the hard axis. Let Hr be the magnitude of 
the rotating field. Let M u be the component 


of the easy-axis (x-axis) irreversible magnetization and M ty be 
the component of the hard-axis (jy-axis) irreversible 
magnetization during rotation of the applied field Then 


M s$x er f 


Hrx-Hcx 


M iy = Ms S y erf 


H Ry ~ H Cy 


where M s is the saturation magnetization, S x and S y are the 
easy and hard axis squareness respectively, H Rx and H Ry are 
the easy and hard axis average critical fields respectively, and 
<j x and o y are the easy and hard axis standard deviation of the 
interaction and critical fields respectively. Define the 
reversible magnetization function, 

G=a + f{H R )-a_f{-H R ) (22) 

The coefficients, a+ and a . in (22) are derived for the easy and 
hard axes using the irreversible magnetization of (20) and (21), 
1 + M ix 1-Mfc 

a x+ = and a x _ = (23) 


a y+ 


and ciy_ = - 


1-M*. 


The reversible function, / in (22), may be any convenient 
function that goes through zero when the applied field equals 
zero and is finite for large values of applied field. Consider 
the function 

/ (H )= 1 -exp(-£l? ) for 7/ > 0 (25) 

and 

f(H)=H for II < 0. (26) 

where £ is a model parameter that needs to be identified. The 
reversible magnetization is 

M R =G(l-S)H R M s 1 Hr , (27) 

where 1 jj r is a unit vector in the direction of H R . For the easy 
and hard axes (x and y respectively), 


=G x (l-S x >f s 


Mjfcrd 



Frg. 1 1 . Vector sum of magnetization contributed by oscillating 
magnetization (Mosc) and magnetization contributed by rotating 
magnetization (M ra ). 


and 

M R y =Gy^~Sy]M S . (29) 

Adding (20) and (23) 

M x = M S S x erf f Hrx ~ | + G x (l - S x )M S ■ (30) 

l ° x 

Adding (21) and (24) 

M y =M s SytTi — ^ +Gj,^-Sj)w5. (31) 

<Jy 

v y 

Then, the ellipticity correction factor, equals (30) divided 
by (31) 
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The rotational magnetization is computed from 



and when added to M osc computes the resultant magnetization, 
M. Multiplying (30) by C 0 and (31) by C r enables the 
computation of 0 H and 0 M , the angle of applied field and the 
angle of magnetization respectively for the rotation of the 
applied field. Figure 12 shows the rotation coefficient, Q , for 
positive applied fields. 

i| , , , , — , 

0.9 • / 

0.8 ■ / 

0.7 • / 

o- / 

1 06 / 
o / 

IE / 

§ 0.5 - / 

o / 

■I 0.4 - / 

o / 

a 0.3 • / 

0.2 • / 

0.1 ■ / 

°0 0.05 0 1 0.15 0.2 D25 0.3 0.35 

Applied Field PqH, T 

Hg. 12 Computed rotation coefficient, Cr , for positive applied fields. 


material is different for each different applied field. The 
magnetization trajectories in Fig. 13 are for applied fields with 
a fixed magnitude and for rotations from 0° to 360°. The 
difference between the starting and ending magnetization 
magnitude is a measure of the accommodation. The difference 
between the starting and ending magnetization angle is a 
measure of the lag. As the applied field increases, the shape of 
the trajectories changes exponentially from ellipses to circles 
as shown in Fig. 13. 

The details of the lag angle and the accommodation are 
shown in Fig. 14 for a measured magnetization MPT trajectory 
for a clockwise rotating applied field of 0.2 T. A 



27D 


Fig. 13. Change in magnetization for applied fields of po H = 0.2 T, 
0.3 T, and 1.0 T 


vm. Rotational Accommodation 

Accommodation, based upon the magnetic interaction 
between hysterons, causes the drifting of minor loops with the 
reversal of applied fields. This drifting has an important 
impact, because it affects the repeatability of the recording 
process and causes performance changes during the operation 
of motors and transformers. To understand accommodation, 
one must understand to where the magnetization is drifting. A 
magnetic calculation in a complex process will give a wrong 
answer with the classic Preisach model that does not have 
accommodation. Accommodation, which is rate independent, 
appears to be similar to aftereffect, which is rate dependent. 
Accommodation and aftereffect are caused by different 
processes. 

Vector measurements made on MP tape have displayed 
vector accommodation [11]. For these measurements, several 
different applied fields were rotated 360° after initialization of 
a negative saturation field. The jc- and y- components of the 
magnetization were measured. The vector sum of these 
measurements, which we call magnetization trajectories, 
generate elliptic loci [11]. A magnetization trajectory of a 
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more precise measure of accommodation is computing the 
difference between starting and ending of the trajectory 
magnitude after the ending magnetization is extrapolated to the 
same angle that the magnetization started (0° and 360°). 


Accommodation 



Fig. 14. Measured MPT magnetization trajectory for a 
clockwise rotating applied field of 0.2 T showing a measure of 
magnetization accommodation and lag. 


9(1 02 



Fig. 15. Measured and Computed Ellipsoidal Trajectory at po H = 0.2 T 


Figure 15 shows the measured MPT magnetization 
trajectory and a computed ellipse with and exponential decay 
for a rotating applied field of 0.2 T. The computed trajectory 
shows that the measured trajectory is ellipsoidal and that the 
accommodation is an exponential decay of this trajectory. The 
decay tends towards zero as the applied field increases to large 
values as shown in Fig. 13. 


IX. Conclusions 

A simpler vector model is developed that does not require 
the rotational correction for computation speed. The 
elimination of the rotational correction makes it possible to 
implement the differential method, a very effective way to 


compute the magnetization. The simulation results show that 
the presented model represents isotropic media accurately. 
Further refinement of this vector model in terms of speed and 
accuracy is a topic of further research. 

The orientation correction for isotropic material was 
extended for anisotropic material. This extension consists of 
identifying rotational and oscillatory magnetization by 
determining the percentage of the hard axis Preisach function 
swept by the applied field. The extension corrects the 
magnetization to produce the anisotropic double fold 
periodicity. The amplitude of the double fold periodic 
magnetization tends towards zero as the applied field 
increases. 

Vector measurements made on magnetic particle tape, and 
anisotropic medium, display vector accommodation, which is 
especially evident in polar format. For small applied fields, 
the magnetization trajectories form ellipses that exponentially 
form circles as the applied field increases to large values. A 
measure of the accommodation is the difference between the 
start and the end of a magnetization trajectory. The 
accommodation may be represented as an exponential decay of 
an ellipsoidal magnetization trajectory. Incorporation of 
vector accommodation into a Preisach vector model is the 
subject of future research. 
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TABLE L SVPM State Function Values in Two Dimensions, Q, 


1 H 


EBB! 

h x >«x 



0 

jMMDBMNHI 


m 

i 

EBB 

B 


§ 

v y <h y <u y 

-1 

no change 

l 

hy >Uy 

MBm 

0 



^iH! 

■ 

EBB! 

B 

Si 

1 


TABLE II RVPM State Function Values in Two Dimensions, Q, 



v x > h x 

v x <h z <u x 

h x >u x 

m 

h x ~v x 

0 


■J( h x~Vxf +^y-V y f 

^(h X -U x f +(fly-Vy} 

V y < hy < Uy 

-1 

No change 

1 

ky >Uy 

h x~ v x 

0 

h x -“x 

^ X -V X f+(fly~Vy} 

■j(hx~U Z f +(fly-Vyf 


table m. State Function Values in three Dimensions, Qj 


Number of 
violations 

Violations 

States 

0 

V • </l <u . 

J V “ J 

holds for j = x,y and z 

no change 

1 

hj >UjOT hj < v j 

Qj =lor Qj =-l and g, =0 ,i* j 

2 

Any two combinations of 
violations in u or v where 
the thresholds violated 
are called tj mdt k . 

_ hj - tj h k —t. 


and Q t = 0 , i * j, k 

3 

Any three combinations 
of violations in u or v 
where the thresholds 
violated are called 
ti,tj andf t . 

n - h ‘~ t( 

" Vfc-cMVLMA-rJ 

J( h i -*i f + {hj -tj J + {h k -t k f 

n - hk ~ tk 

^{h, ~ t, ) 2 + {hj -tjj + (fc t - t k f 







































